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System of 1/2 spin particles is observed repeatedly using 
Stern-Gerlach apparatuses with rotated orientations. Synthe- 
sis of such non-commuting observables is analyzed using max- 
imum likelihood estimation as an example of quantum state 
reconstruction. Repeated incompatible observations repre- 
sent a new generalized measurement. This idealized scheme 
will serve for analysis of future experiments in neutron and 
quantum optics. 



I. INTRODUCTION 

Quantum mechanics of 1/2 spin particles serves often 
as an illustrating example for many quantum considera- 
tions in standard textbooks of quantum theory |^] . The 
importance of this model is enhanced by the fact that 
such states represent the smallest amount of quantum in- 
formation - quantum bits (q-bits). Besides theoretically 
valuable "Gedanken" experiments, the spin 1/2 particles 
such as electrons, neutrons or polarization states of light 
quanta are convenient for feasible experiments in mat- 
ter wave and quantum optics. They play crucial role 
in many sophisticated schemes involving entanglement, 
Bell state analysis or teleportation. Several approaches 
for measurement and estimation of spin states have been 
considered recently |^|^. In this Brief Report, the max- 
imum likelihood (MaxLik) estimation of 1/2 spin state 
will be formulated as an illustrating example of more 
general treatment The given formulation shows 

a tight relation between quantum and statistical theo- 
ries. Synthesis of many independent and nonequivalent 
ideal detections of Stern-Gerlach (SG) type will be in- 
terpreted as a new generalized measurement, output of 
which the quantum state is. This consideration will be 
used as a basic tool for further investigation in depolar- 
ization measurements, neutron and light interferometry 
and quantum state reconstruction. 

Basic properties of spin 1/2 quantum systems will be 
briefly reviewed. A pure state (projector) may be repre- 
sented by the expression 



(1) 



where a = (ai,a2,a3) is the three-dimensional normal- 
ized vector, ai,i — 1,2,3 represent the Pauli matrices and 



the summation convention for repeated indices is used. 
Since 



aia^ — 6ij + iei-jkCTk, 



the scalar product of two projectors is given as 

1 



\{a\b)\' 



-(1 + aib,). 



General state described by a density matrix may be pa- 
rameterized by 



P = P+|a)(a| +p-|-a)(-a| 
- + -aiai{p+ -P-), 



(2) 
(3) 



where p+ + — 1 and the states |a),|— a) denote a 
general orthogonal basis. Spin state may be alternatively 
described by a polarization vector 



(4) 



where the brackets () denote an expectation value. Hence 
the polarization completely determines the state of 
quantum system. Degree of polarization may be intro- 
duced as 



IrP < 1 



and |rp = for completely unpolarized (mixed) state 
and |rp = 1 for completely polarized (pure) states. 

The polarization or spin may be measured by project- 
ing the state into the given directions of SG apparatus 
±a. Closure relation and operator representation of such 
a device simply read 



A 



|a)(a| + |-a)(-a| = l, 
:i[|a)(a|-|-a)(-a|]. 



(5) 
(6) 



Assuming for the sake of simplicity always the same to- 
tal number of particles N, the number of particles with 
the spin "up" and "down" estimates projections of the 
polarization vector according to the relations 



n± = Np{±a.) = ^N{1 ± ra). 



(7) 



Since this may be done for three orthogonal directions 
of coordinate axes Xj, i = 1, 2, 3, the polarization may be 
found by eliminating the total number of particles N 
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H+ 



(8) 



Each polarization component is determined separately. 
This represent a correct solution, provided that resulting 
polarization is inside the Poincare sphere |rp < 1 only. 
This example has been used for motivation of general 
analysis of quantum state reconstruction in Ref. ||^. As 
demonstrated, the "states" outside the Poincare sphere 
violate the positive semidefiniteness of quantum states 
yielding improper quantum description of noises. Similar 
problems appear in the case when more than three pro- 
jections are used. Some results of SG projections might 
appear as incompatible among themselves due to the 
fluctuations and noises involved. Various SG measure- 
ments are not equivalent, since they are observing vari- 
ous "faces" of the spin system. Such measurements, even 
when done with equal number of particles, determine var- 
ious projection with different errors. Detected data ni ± 
collected from SG observations in M various directions 
±a', i — 1,2,... Ad sample various binomial distribu- 
tions. Significantly, the detected data ni^± fluctuate with 
the root -mean square errors given by y/ N{1 — (raJ)^)/2 
depending on the deviations between projections and the 
true (but unknown!) direction of the spin r. Various pro- 
jections cannot be therefore trusted with the same degree 
of credibility, since they are affected by different errors. 
This is manifested in quantum theory, since various SG 
measurements are incompatible and corresponding oper- 
ators (|l]) do not commute for different orientations a^ . 
Such data cannot be obtained in the same measurement, 
but may be collected by repeating. An optimal proce- 
dure must predict an unknown state and simultaneously 
take into account data fluctuations. This indicates the 
nonlinearities of an algorithm. MaxLik estimation does 
this job and fits the data to a quantum state. Besides 
this, it is the only procedure, which provides the same 
structure as generalized measurement [^j. Henceforth, 
synthesis of incompatible measurements may be inter- 
preted as a novel measurement of quantum state. This 
will be demonstrated in the following section. 



II. SPIN ESTIMATION 

Provided that source supplies 1/2 spin particles pre- 
pared in the same mixed state, an ideal lossless SG mea- 
surement performed repeatedly on the system of TV par- 
ticles will be assumed. The setting of SG apparatus may 
change. Provided that detection has been done with M 
different settings, N * M particles have been used alto- 
gether and an unknown quantum state should be found. 
The results of the measurement may be characterized by 
settings of the SG apparatus ±a-' and by the relative 
frequencies of the outcomes 1/2(1 ± Xj) — nj^±/N. The 
question is what state(s) fit(s) the data in optimal way. 
One might be tempted to sample and invert the proba- 
bility, predicted by quantum theory, as it is done in the 



case of equation (||). Because each SG detection is repre- 
sented by a complete measurement, the sum the relations 
(H) for each setting of SG apparatus j reads 



M 



However, the expected relations 

Tr{p— |±aJ)(±aJ'|} = — (1±X,). 



(9) 



(10) 



cannot be fulfilled, in general, since the system is over- 
completed and data are fluctuating. Hence, the proba- 
bilities cannot be mapped so straightforwardly with the 
relative frequencies of outcomes. 

MaxLik principle provides a tool, how to treat this 
problem. The most probable state consistent with the 
data should be found. As the measure of probability, the 
likelihood functional corresponding to the product of all 
the probabilities for all detected data may be constructed 



W(l+X,)/2 



-ar' |/5|— a^ 



N{l-X,)/2 



(11) 



Extremal states of likelihood functional satisfy the non- 
linear operator equation 0| 



2M ^ 
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(1+^.) 



|a^")(a^l 
{ao\p\ai) 



(1 - X, 



(-aJ|p|-a^ 



P = P 



(12) 



Quantum state may be represented by polarization. Us- 
ing the relation (^ , multiplying both the sides by Ufc and 
performing the trace, the equation reads 



i?(r)r -I- K(r) + iK(r) x r = r, 
where the functions are defined as 



(13) 



^1 2M ^ V 1 
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1 - X, 



a>v 1 a^ r 



X, 



l-X, 



a?v 1 aJ r 



Because real and imaginary parts are dependent, the real 
part of this equation represents sufficient and necessary 
conditions. The final equation for polarization reads 



i?(r)r + K(r) = r. 



(14) 



The relation ([l4| ) may be advantageously used for itera- 
tion of the solution. Starting from the centre of Poincare 
sphere r = 0, the left hand side of eq. ( [1^ ) provides the 
first correction, which may be used for subsequent itera- 
tion, etc.. This procedure provides a quick algorithm for 
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MaxLik fitting of an unknown quantum state inside the 
Poincare sphere. 

An equivalent result may be derived, provided that 
likelihood function is parameterized directly using the 
polarization. The relevant part of the likelihood function 
corresponding to the observation of particular data reads 

/:((r)) = 11(1 + raJ)^(i+^.)/2(i _ ra^-)^(i-^.)/2. (15) 
j 

The vector r parameterizes an unknown polarization in- 
side the Poincare sphere and the products runs over all 
M directions. The standard statistical approach using 
MaxLik ^ In £ provides the vector equation for extremal 
polarization 

yZl_^a^^O. (16) 

j 

The equation (|^ is equivalent to the equation (|l6|). In- 
deed, the equation (|l^) is nothing else as K(r) = 0, im- 
plying the relation i?(r) = 1. On the other hand, the 
equation (|lj) may be rewritten to the form of equation 
(|^) as well. 

Results of numerical simulation are shown in the Fig. 
|l|. Stern-Gerlach detection is simulated here for projec- 
tion of an unknown state (north pole on Poincare sphere) 
in five various directions. Each measurement is done with 
20 impinging particles, which are registered either with 
the spin up (upper left panel) either down (lower left 
panel) . Both the left panels show typical values for a sin- 
gle experiment. For each position of the projector there 
are three bars in the upper and lower left panels. The first 
bars (black) show the true value of the probability. The 
second bars (gray) show the counted statistics fluctuating 
around the true value of probability. The third bars (hol- 
low) show results of the reconstruction - the statistics of 
reconstructed state corresponding to the given projector. 
Notice here, that upper and lower panels are comple- 
mentary and sum of corresponding probabilities is one. 
The right panels visualize results of 10 times repeated ex- 
periment on the Poincare sphere. Symbols of diamonds 
denote the positions of five projectors. Orthogonal pro- 
jectors in opposite directions are not depicted. The stars 
show the position of reconstructed states. The true state 
corresponds to the north pole. Lower right panel show 
the upper view. 

Quantum formulation posses a nontrivial interpreta- 
tion, which can hardly be recognized in the equation for 
polarization (p^). Because the above scheme determines 
a quantum state, it must exist a generalized measure- 
ment described by a probability operator measure (POM) 
result of which the quantum state is. Really, such a 
probability operator measure can be easily find by proper 
renormalization of original SG measurement Let us 
define the POM as rcnormalized SG projectors 

\±a^){±a^R = ---i^^^_|±a^')(±a''l, (17) 
' ' 2M(±a-J|pe|±a-') ' ^ ^ 



for each index j. The closure relation then reads 

M 

^|a^)(a^U + |-a^")(-a^U-i, (18) 
i 

and the rcnormalized POM fulfills the conditions 

Tr{pe.*|±a^)(±a^k} = ^(1 ±^.)- (19) 

Here pe denotes the extremal state - a solution of the 
equation ([l^). Indeed, the relation ( |l8|) coincides with 
the equation for extremal states (p^, whereas the con- 
dition for expectation values (|l9|) is fulfilled as an iden- 
tity. The reconstruction is done on the subspace, where 
the rcnormalized POM reproduces the identity operator. 
Particularly, this means that the identity operator on the 
right hand side of ( |l8| ) is spanned by the one dimensional 
subspace only (i.e by a single ray), provided that the ex- 
tremal state Pe is a pure state. For a general extremal 
density matrix, the reconstruction is accomplished in the 
whole two dimensional Hilbert space. The distinction 
between relations (I), (0) and (|l|), characterize 
the subtle point of quantum state reconstruction. Max- 
Lik solution may be also interpreted in the language of 
probabilities. The detected data ni_± samples different 
binomial probability distributions for i — 1, . . . , M. Max- 
Lik estimation finds a common multinomial distribution, 
sampling of which the data seems to be with the highest 
likelihood. 

The method developed here may be compared with 
the existing approaches. Jaynes maximum entropy prin- 
ciple (MaxEnt) ||l^ has been applied to the estimation 
of spin 1/2/ states in Refs. |^,^. In general, these meth- 
ods are not equivalent. The MaxLik solution seeks for the 
most likely solution consistent with the data, whereas the 
MaxEnt searches for the worst solution still consistent 
with the data. This may be interpreted as different prior 
information in the maximum probability principle Jll| . 
But this is not the only difference. External conditions 
of both the approaches differ substantially. MaxLik has 
been applied to the measurements with many projectors. 
However, the same conditions cannot be used in MaxEnt 
principle. Since there are only three free parameters nec- 
essary for determination of an unknown spin state, con- 
ditions ( p^ ) cannot be fulfilled, in general. MaxEnt prin- 
ciple cannot be used provided that there are more than 
three independent conditions put on the density matrix 
of 1/2 spin system. 

In the papers Refs. |l^ an optimal strategy for mea- 
suring an unknown two-state system prepared in a mixed 
state is investigated. As a result, an optimal POM may 
be predicted. On the other hand, not the measurement 
but the mathematical treatment is optimized in this pa- 
per. This seems to be reasonable from the experimen- 
talist's point of view since it is questionable how to do 
a general measurement described by a POM. As demon- 
strated, for given measurement MaxLik estimation pro- 
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vides an optimal treatment, since it reproduces a quan- 
tum measurement. As argued in Ref. |l3[ both the treat- 
ments provide comparable results. 



III. SUMMARY 

Synthesis of incompatible observations represented by 
various settings of SG apparatus has been evaluated us- 
ing MaxLik estimation. As the result it defines a gen- 
eralized measurement of a quantum state. Synthesis of 
various incompatible observations is again a generalized 
measurement. Developed formalism will be applied in the 
future for investigation of various problems such as: spin 
estimation of neutrons in the depolarization experiments, 
estimation of quantum state inside the interferometer or 
analysis of entanglement. 
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FIG. 1. Results and interpretation of spin reconstruction 
for numerical simulation of Stern-Gerlach detection with ap- 
paratuses in five various settings 
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